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Abstract: The scaling form of the free-energy near a critical point allows for the definition of various 
universal ratios of thermodynamical amplitudes. Together with the critical exponents they character- 
ize the universality classes and may be useful experimental quantities. We show how these universal 
quantities can be computed for a particular class of universality by using several Quantum Field Theory 
methods 



1. Introduction 

One of the most important successes of Quantum 
Field Theory (QFT) in recent years consists of 
the quantitative analysis of the universality classes 
of two-dimensional statistical mechanics models 
near their second order phase transition points. 
Right at criticality, where the correlation length 
£ diverges, Conformal Field Theory (CFT) allows 
the exact determination of the spectrum of anoma- 
lous dimensions, the structure constants of the 
OPE algebra, the multi-point correlation func- 
tions etc. [Q. However, the CFT data do not ex- 
haust all physical information relative to the phase 
transitions. From a theoretical point of view, a 
perturbation of the conformal action is required 
in order to investigate the space of the coupling 
constants and its topology consisting of the loca- 
tion of the fixed points and the Renormalization 
Group flows which connect them. This rich and 
interesting subject has been deeply investigated 
by several authors in the last ten years (see, for 
instance, || and references therein for a review). 
From a practical point of view, one has to take 
into account that in a real sample, conformal in- 
variance may be broken by the presence of impu- 
rities or, simply, by an imperfect fine-tuning of 
the experimental knobs. Hence, to control these 
effects it becomes important to study the dynam- 
ics of the systems slightly away from criticality, in 
particular their responses to external fields. The 
most ambitious goal would be the determination 
of both the scaling function which fixes the equa- 



tion of state and the off-critical correlators of the 
various order parameters. Although the exact de- 
termination of the equation of state of a given uni- 
versality class may often be a difficult task, the 
scaling property alone of the free-energy is never- 
theless sufficient to extract numerous predictions 
on universal combinations of critical amplitudes. 
As will be discussed below, these universal combi- 
nations are pure numbers which can be extremely 
useful for the experimental identification of the 
universality classes. 

2. Universal Amplitude Ratios 

Consider a statistical model with n relevant fields 
ipi(x) at criticality. Near the critical point, its 
action can be parameterised as 

A = AcFT + gi J ¥i{x)d D x . (2.1) 

The scaling property of the order parameters is 
encoded into the asymptotic form of their two- 
point functions {(fii(x)ipi(0)) ~ , h% . for |x| — > 
so that Ai is identified with the conformal dimen- 
sion of the fields. Correspondingly the conjugate 
coupling constants <?.; behave as gi ~ \D-2&i ^ 
where A is a mass scale. Therefore, away from 
criticality there will be generally a finite correla- 
tion length £ which in the thermodynamical limit 
scales as £ ~ a (Kigi)~ -°- 2A ; , where a ~ A^ 1 
may be regarded as a microscopic length scale. 
The terms Ki are non-universal metric factors 
which depend on the unit chosen for measuring 
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the external sources gi, alias on the particular re- 
alization selected for representing the universality 
class. In the presence of several deformations of 
the conformal action, the most general expression 
for the scaling form of the correlation length may 
be written as 



£ = & = a{Kigi)~ 



(2.2) 

up = D _ 2A are the so-called crossover 
exponents whereas Li are universal homogeneous 



where 



Z5-2A, 



K i9i 



- . There are 



scaling functions of the ratios 
of course several (but equivalent) ways of writing 
these scaling forms, depending on which coupling 
constant is selected as a prefactor. In the limit 
where gi — > (I =fc i) but gi ^ 0, equation ( |2.2| ) 
becomes 



(2.3) 
g n \. This 



Consider now the free-energy f\g\ 
is a dimensionless quantity defined by 

e -f{9i,~,9n) _ J x>0 e -[^FT+X;r=i9> J Vi(^)<i D ^ 

(2.4) 

Assuming the validity of the hyperscaling hypoth- 
esis, in the thermodynamical limit its singular part 
(per unit of volume) will be proportional to the D 
power of the correlation length. Let us denote the 
singular part of the free-energy for unit volume by 
f[g±, . . . , g n \. Depending on which scaling form is 
adopted for the correlation length, we have cor- 
respondingly several (but equivalent) ways of pa- 
rameterizing this quantity 



fi[gi,...,g„] = (K igi )T^i Ti 



{K i9i )o 



(2.5) 

The functions J~i are universal homogeneous scal- 
ing functions of the ratios - — . 

Let us consider now the definition of the ther- 
modynamical quantities related to the various deriva- 
tives of the free-energy. We will adopt the nota- 
tion (...}{ to denote expectation values computed 
in the off-critical theory obtained by keeping (at 
the end) only the coupling constant gi different 
from zero. The first quantities to consider are the 
vacuum expectation values (VEV) of the fields ipj 
which can be parameterized as 

dfi 



9 9j 



= B j% g 



2A j 

D-2A; 



(2.6) 



with 



B 4 



2A 3 



><■ ~i\<*; " ■ (2-7) 

The above relations can be equivalently expressed 
as 

9i = Aj 



with 



Di 



1 



(2.8) 
(2.9) 



The generalized susceptibilities of the model are 
defined by 



r 



jk 



_d_ 

dgk 



d 2 f t 



dgkdgj 



(2.10) 



Sl=0 



They are obviously symmetrical in the two lower 
indices. By extracting their dependence on the 
coupling constant gi, they can be expressed as 



— r}fc g. 



2A j+ 2A fc--D 
D-2A; 



with 



2A i +2A l -D 



r 



jk 



(2-11) 



(2-12) 



The various quantities obtained by taking the 
derivatives of the free-energy obviously contain 
metric factors (the quantities Ki) which make their 
values not universal. However, it is always pos- 
sible to consider special combinations thereof in 
such a way to cancel out all metric factors. Here 
we propose the consideration of the following uni- 
versal ratios 



{Rc 



'jk 



T^2 pi 

ii jk 

BjiBki 



(R x )) = TjjDjjBji 

4 = (4) 1/D ^ 

4A-+2A;-2D 2 



(R A )) = r) j D ii B-, A 



(Q 



' 30 it 

r 



2)jk - 



j j 



r 



D-4A, 



(2.13) 

(2.14) 
(2.15) 

(2.16) 
(2.17) 



9l=0 



33 \~=3 / 

From their definition, these quantities are pure 
numbers attached to the universality classes and 
can be used to characterize them. In fact, the am- 
plitude ratios are numbers which typically present 
significant variations between different classes of 
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universality, whereas the critical exponents usu- 
ally assume small values which only vary by a 
small percent by changing the universality classes. 
Hence the universal ratios may be ideal marks of 
the critical scaling regime It is also worth 
emphasizing that, from an experimental point of 
view, it should be simpler to measure universal 
amplitude ratios rather than critical exponents: to 
determine the former quantities one needs to per- 
form several measurements at a single, fixed value 
of the coupling which drives the system away from 
criticality whereas to determine the latter, one 
needs to make measurements over several decades 
along the axes of the off-critical couplings. More- 
over, although not all of them are independent, 
the universal ratios are a larger set of numbers 
than the critical exponents and therefore permit 
a more precise determination of the class of uni- 
versality. 

3. Tricritical Ising Model 

The universal ratios of the two-dimensional Tri- 
critical Ising Model (TIM) were computed in refs. || 
0]. We refer the reader to these articles for a de- 
tailed discussion of their evaluation. There are 
several reasons for considering the class of uni- 
versality associated to this model. First of all, 
from the experimental point of view, a number of 
physical systems exhibit a tricritical Ising behav- 
ior, among them fluid mixtures, metamagnets or 
order-disorder transitions in absorbed systems 
(for a review on the theory of tricritical points, see 
Q). Secondly, from a theoretical point of view, 
the TIM is still sufficiently simple to be solved 
but at the same time it presents an extremely 
rich and fascinating structure of excitations away 
from criticality. In a Landau-Ginzburg approach, 
the TIM is associated to a 3> 6 -theory near its 
tricritical point |^], although this approach is of- 
ten too elementary for the understanding of some 
of its remarkable symmetries. Depending on the 
direction in the phase space in which the sys- 
tem is moved away from criticality, one can ob- 
serve, for instance, a behavior ruled by the excep- 
tional root system E 7 jl^, |ll| or by supersymme- 
try [jl^L [l3[ [l4|, |l5| (in its exact or broken phase 
realization) or by an asymmetrical pair of kinks 
iPI [It], In addition, the description of its 



low-temperature phase is easily obtained from the 
one of its high-temperature phase because of the 
self-duality of the model. 

A convenient way to determine the universal 
ratios of the TIM consists of adopting a Quantum 
Field Theory (QFT) approach, as it was done in 
the original references || |(| . It would be obviously 
interesting to compute and to compare them by 
using a lattice formulation as, for instance, the one 
studied in 1 1 9 . In a QFT approach, one takes ini- 



tially advantage of the exact solution of the model 
at criticality. The bidimensional TIM is described 
by the second representative of the unitary series 
of minimal models of CFT |) : its central charge 
is equal to c = and the exact conformal weights 
of the scaling fields are given by 



Aj.fc = 



(5Z -4fc) 2 



80 



1 < I < 3 
1< k < 4 



(3.1) 



There are six primary scalar fields (j> A A , which 
close an algebra under the Operator Product Ex- 
pansion 

4>i(zi,zi) (pj(z 2 ,Z2) ~ (3.2) 

~ £ C ijk \ Zl -Z2 |-2(A i+ A,-A fc ) fafaz*) . 
k 

The skeleton form of this OPE algebra and the 
relative structure constants of the Fusion Rules 
of the TIM were given in ref. The six pri- 

mary fields can be identified with the normal or- 
dered composite LG fields 0. With respect to 
their properties under the Zi spin-reversal trans- 
formation Q : $ — » — $ we have: 

1. two odd fields: the leading magnetization 
operator a — 6 3 3 = $ and the sublcading 

80 I 80 

magnetization operator a' = 6 i i =: ^ : 

° 16 ' 16 

2. four even fields: the identity operator 1 = 
</>o,o> the leading energy density e = </>_l a. ='■ 
<I> 2 :, the subleading energy density t = </>e. t e_ 

:, which in metamagnets assumes the 
meaning of the density of the annealed va- 
cancies, and the field e" = 63 3 . The OPE 

2 ■ 2 

of the even fields form a subalgebra of the 
Fusion Rules. 

In the TIM there is another Z2 transformation - 
the Kramers- Wannier duality D - under which 
the fields transform as follows: 
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• the order magnetization operators are mapped 
onto their corresponding disorder operators 



D-WD 



(3.3) 



• the even fields are mapped onto themselves, 
D~ l eD = -e , D~ 1 tD = t , (3.4) 

Given the ^-spin odd parity of the two magnetic 
fields, the off-critical theories obtained by their 
perturbation are independent of the sign of their 
associate couplings. The odd parity of the energy 
field under the other Zi symmetry is responsable 
for the duality mapping between the high and low 
temperature phases. On the other hand, two dif- 
ferent physical situations arise depending on the 
sign of the coupling relative to the perturbation of 
the vacancy density field t(x). It is important to 
notice that, excluding the leading magnetic per- 
turbation, all the others give rise to integrable 
QFT away from criticality. To simplify the formu- 
lae below, it is convenient to adopt the compact 
notation ipi (i — 1,2,3,4) to denote collectively 
all these fields, so that tpi = a , P2 = £, <Pz = o' 
and tp^ = t. Let us discuss now some aspects of 
QFT which make possible the determination of 
the universal ratios and the strategy we have used 
to achieve these results. 

4. Quantum Field Theory Approach 

Each coupling constant (i = 1, ... ,4) relative 
to the relevant operator ipi(x) of the TIM is a di- 
mensional quantity which can be related to the 
lowest mass-gap m, = C,^ 1 of the off-critical the- 
ory according to the formula 



mi = Ci g- 



(4.1) 



When the QFT associated to the action (2.1) is in- 
tegrable, the pure number Ci can be exactly deter- 
mined by means of the Thermodynamical Bethe 
Ansatz [|o[ |l| . When the theory is not integrable 
(this is the case for the magnetic deformation of 
the TIM), the constant Ci can be nevertheless de- 
termined by a numerical method, based on the so- 
called Truncated Conformal Space Approach p2| . 
In conclusion, for all individual deformations of 
the TIM we are able to completely set the rela- 
tionship which links the coupling constant to the 



mass-gap of the theory and therefore to switch 
freely between these two variables. 

Another set of quantities which can be fixed 
by QFT are the matrix elements of the order pa- 
rameters, the simplest ones being the vacuum ex- 
pectation values (VEV). In this case we have 



(tpj)i = B 



31 ; 



(4.2) 



When the theory is integrable, the constant Bji 
can be fixed exactly, thanks to the results of a 
remarkable series of papers [^3[ ^4|. When it is 
not integrable, the constant Bji can be neverthe- 
less estimated by means of a numerical approach, 
as firstly shown in |25| . Hence, also in this case, 
we are able to determine completely these quan- 
tities. Moreover, as shown in a generaliza- 
tion of the numerical approach of ref. often 
leads to a reasonable estimation of the matrix ele- 
ments of the order parameters between the vac- 
uum states and some of the excited states, as 
for instance (0|<^j|j4.fc)i where Ak is a one-particle 
state of mass Mfe. These quantities turn out to be 
useful for obtaining sensible approximation of the 
large-distance behavior of several correlators. 

Another useful piece of information on the off- 
critical dynamics can be obtained by exploiting 
the properties of the stress-energy tensor T^ v {x). 
In the presence of the perturbing field pi , the trace 
of the stress-energy tensor is different from zero 
and can be expressed as 



9(x) = 27r ffi (2-2Ai)<pi 



(4.3) 



The trace of the stress-energy tensor enters a use- 
ful sum rule - called the A-theorem sum rule |3l| 
- which reads 



A ; 



1 



An{ipj 



d 2 x(Q(x)p J (0)) c l > (4-4) 



i.e. it relates the conformal dimension Aj of the 
field <pj to its VEV and to the integral of its con- 
nected off-critical correlator with <d(x). It is easy 
to see that the above formula simply expresses 
the content of the fluctuation-dissipation theo- 
rem and when the above integral diverges, so does 
the VEV in the denominator, in such a way that 
eq. (4.4) always keeps its validity ||. 

Basic quantities in the universal ratios are the 
generalized susceptibilities r* fe which, by using the 
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fluctuation-dissipation theorem, can be expressed 
as 

Y) k = J d 2 x(<p 3 (x)<p k (0)y c . (4.5) 
By extracting their dependence on the coupling 



constant gi, we have f = r*. fc g i 



l-A; 



with 



(4.6) 

Some of the above susceptibilities can be deter- 
mined exactly, such as the components T\ k , whose 
values are provided by the A-theorem sum rule 



Aft 



1 



-Bi 



(4.7) 



In all other cases, when an exact formula is not 
available, the strategy to evaluate the generalized 
susceptibilities relies on two different representa- 
tions of the correlators. These representations 
have the advantage to converge very fast in two 
distinct regions. 

The first representation is based on Conformal 
Perturbation Theory (with the employment of the 
non-analytic expression of the VEV) |52|. In this 
approach the two-point correlators are expressed 
as 



\CUg;x){A p {Q)) (4. 



where the off-critical structure constants Cj fe (g; x) 
admit the expansion 

oo 

Cf j (g;x)=r 2 ^- A ^ £ C^\gr^*Y , 

71=0 

(4.9) 

(r =| x |) and can be computed perturbatively in 
g. Their first order contribution is given by |32| 



(4.10) 

where the prime indicates a suitable infrared (large 
distance) regularization of the integral. This rep- 
resentation allows a very efficient estimation of the 
correlation function in its short distance regime 
r < f. 

The second representation is based on the Form 
Factors and allows an efficient control of its large 
distance behavior, i.e. when r»(. In this second 
representation, one makes use of the knowledge of 



the off-critical mass spectrum of the theory to ex- 
press the correlators as 



(<Pi(x)<pj(0)) = J2g n (r 



(4.11) 



n=0 



where 
9n(r) = 



d0\ dO n _ r V™ rrifccosh* 



2vr 2vr 



H>e 3 ...>e n 

x(0\ip i (0)\...A^{9 n ))(...A an {6 n )\(p j (fl)\0) . 

\A ai (0i) . . . A an (#„)) are the multi-particle states 
relative to the excitations of mass m k , with rela- 
tivistic dispersion relations given by E = m k cosh 8 , 
p = TOfcsinh#, where 9 is the rapidity variable. 



The spectral representation (4.12) is obviously an 



mi 



expansion in the parameter e~ f , where £~ 
is the lowest mass-gap. 

Basic quantities of the large distance approach 
are the Form Factors (FF), i.e. the matrix ele- 
ments of the operators (fi on the physical asymp- 
totic states |26l Wft, 



JTVi 
ax, 







01^(0)14,, A an (9 n )) 
(4.12) 

It is worth emphasizing that the above quanti- 
ties are unaffected by renormalization effects since 
physical excitations are employed in their defini- 
tions. For scalar operators, relativistic invariance 
requires that the FF only depend on the rapidity 
differences 9i — 9j . In our calculation we have only 
inserted into the spectral representations the one- 
particle and two-particle FF, computed according 
to the analysis of ref . |p8[ ^9) . 
Both the representations 



and ( 1.11 ) are 

known to converge very fast (see, for instance |3C| , 
|32f and therefore they are efficiently approximated 
by their lowest terms, which therefore can be eval- 
uated with a relatively little analytical effort. These 
considerations obviously lead to the estimation of 
the integral (|4.5|) according to the following steps: 



1. Express the integral in polar coordinates as 



r;. fe = 27r / drr{^(r)<p k (0)y c , (4-13) 
Jo 

and split the radial integral into two pieces 
as 



I = 



r+oo 

/ drr< Vi (r)^(0)>j 
Jo 
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drr(...)l + 



h(R)+I 2 (R) 



+ 00 



dr 



'(•■•>c 

(4.14) 



2. Use the best available short-distance repre- 
sentation of the correlator to evaluate I\ (R) 
as well as the best available estimate of its 
large-distance representation to evaluate h{R)- 

3. Optimize the choice of the parameter R in 
such a way to obtain the best evaluation of 
the whole integral. In practice, this means 
looking at that value of R for which a plateau 
is obtained for the sum of I\ (R) and I2 (R) ■ 
Say in another way, R belongs to that in- 
terval where there is an overlap between the 
short-distance and the long-distance expan- 
sion of the correlator (see, for instance, Fig- 
ure 1, relative to the correlator (ip i(x)<pi(Q))i 
in the high-temperature phase). 



Figure 1: Continuous line = UV approximation, 
dashed line = IR approximation. An overlap of the 
curves is observed around mr ~ 1. 



Gathering all the results relative to the above quan- 
tities, a set of universal ratios for the TIM have 
been obtained. Some of them are exact, like {R c )\ u 
f§A fc , {R c )l k = §A k (k = 1.....4). Those reh 
ative to the low and high temperature phase of the 
model denoted by an upper index =p respectively, 
are in Table 1. An interesting universal ratio is 
provided in this case by the correlation lengths 
^ , measured at the same displacement above and 
below the critical temperature (as extracted from 
the correlation function of the magnetic operator 



using its duality properties) 

f+ /5tt\ 
2- = 2 cos — w 1.28557. 
£- V 18 



(4.15) 



which can be inferred by the exact mass spectrum 
of the model and the parity properties of the ex- 
citations jLOl |l6| |. Other universal ratios are pre- 
sented in Tables 2-5 and a more complete set of 
them can be found in ref . || . 

In view of the predictivity showed by the the- 
oretical approach, it would be interesting to have 
their experimental confirmation. It would be equally 
useful to apply the methods discussed here to other 
models in such a way to bridge a closer contact be- 
tween theoretical and experimental results in two- 
dimensional physics. 
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Table 1: Amplitude ratios R 2 



jk 



i? 2 ! 3.54 


i? 2 3 -2.06 




R-24 = "I 


Rh = 1-30 


i?44 = 1 



Table 2: Universal ratios (R c ))k an d {R( 



)jk ■ 



{R c )\ 2 = 1.05 lO" 2 


(R c )l 3 = 4.85 lO- 2 


{R c )l 4 6.7 lO" 2 


(i? c )33 3.8 1 0- 1 


(i? c )2- = 2.0 10- 3 


(R c )u = -2.34 lO- 2 


{R c )'is 1.79 lO" 2 


(i? c )2 3 " = 3.4 10- 1 



Table 3: Universal ratio (R x )) for i,j — 1,2. 



(R x )\ 3.897 lO- 2 


(i? x )2+ = 0.1111 


(R x )l = 0.116 


{R x f{~ = 0.040 


W = 


{R x )t = 0.1111 



Table 4: Universal ratios i?| and (Ra)) for i,j = 
l,2-,2+. 



Rl 7.557 lO" 2 




R 2+ 1.0784 10- 1 


Rf~ 8.389 10- 1 


(Ra) 2+ 


{R A )l- 3.918 lO- 2 


(R A )l + = 2.958 10- 1 


(Ra)x~ = 8 -260 10- 1 



Table 5: Universal ratios (Q2))k f° r h 3i & = 1, 2 + , 2 



(#2)2+1 


1.260 


(#2)3-1 


= 1.884 


((92)2+2+ 


1.973 


(#2)2+2- 


= 1.320 


Wit 


= 1.56 


(# 2 )?r 


= 0.442 


(Qa)£- 


= 1.70 





